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Abstract 

The field strength correlators at zero temperature are semi-classically evaluated 
fitting the random instanton liquid model to lattice data for quenched SU(3) 
lattice gauge theory. We restrict ourselves to the lowest order in an instanton 
density expansion necessary to explain the difference between D± and T>u. In the 
instanton-instanton and instanton-antiinstanton contributions the Schwinger 
line factors - neglected in a previous analysis - are numerically taken into account 
in a weighted Monte Carlo evaluation. This leads to different estimates for 
instanton size and density. A reasonable description of the correlators within 
the intermediate range from 0.4 fm to 1 fm is obtained. 



I. INTRODUCTION 

Based on the concept of gluon field-strength correlators, a model independent, albeit system- 
atic description of non-perturbative effects in QCD has been worked out. Few correlation 
functions play an important role in the stochastic confinement model and give a detailed de- 
scription of the level splitting of heavy and light QQ bound states M . Relying on the analytical 
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continuation from Euclidean to Minkowski space, they provide a bridge between low and high 
energy QCD phenomenology: they are the basic ingredients for a description of high-energy 
hadron and quark- ( ant i) quark scattering within the stochastic vacuum approach |J. 

By now, numerical results are available from lattice simulations concerning the two-point 
field strength correlators for pure gauge theory with the gauge groups SU (2) |(^-|8j and SU (3) 
[PHl2H over physical distances ranging up to 1 fm. The correlators have also been obtained near 



the deconfinement transition in pure SU(3) gauge theory |]10[ . Somewhat later, this study has 
been extended to full QCD with four flavours of dynamical staggered quarks [fLl] . 

In most of the recent lattice computations of the gluonic field strength correlators the signal 
has been extracted using a cooling procedure which serves the purpose to erase short-range 
fluctuations on the level of a few lattice spacings, thereby redefining (renormalizing) the field 
strength operator on the lattice. Apart from this renormalization, the correlator is expected to 
be not affected by this procedure as long as the diffusive cooling radius r oc \fn remains small 
compared to the distances d of interest (in units of lattice spacing) for an appropriate number 
n of cooling iterations. Presently, this prevents the knowledge of the field strength correlator at 
very short distances below 0.1 fm. Going to these distances requires to do such simulations on 
a few times finer lattices. In the finite temperature case and for QCD with dynamical fermions 
the correlators are presently known only for distances larger than 0.4 fm. 

With cooling like this, eventual semiclassical structures underlying the vacuum of Yang- 
Mills theory or full QCD can hardly be revealed from Monte Carlo lattice configurations. 
Nevertheless, one might try - on the basis of the lattice data for the field strength correlators 
- to justify and to constrain semiclassical models of vacuum structure. By definition, corre- 
lators formed by the field strengths of particular semiclassical configurations do not cover the 
perturbative part of what is measured on the lattice. In order to compare, one has to rely 
completely on the way how the lattice data are split into a (singular at short distances) lowest 
order perturbative part and the non-perturbative signal. This is the way how the renormalon 
ambiguity is presently dealt with (for a discussion see the recent talk by Di Giacomo JT3fl). 



The non-perturbative signal is usually modelled by exponential contributions to the two basic 



structure functions T> and T>\ (see below). These can be replaced by other functional forms 
fl4}1 , e.g. suggested within high energy scattering phenomenology [[15], or provided by some 



specific vacuum model as in the present paper. 

Instantons - localized finite action solutions of the Euclidean Yang Mills field equations - 
are well-known examples of semiclassical configurations [[16] , which have been put into hadronic 
phenomenology in a dilute gas model in Ref. ||17|| . After the IR divergence had been cured 



in a rather ad hoc way taking interactions into account [[Tj|], the instanton liquid model has 
been developed by Shuryak |19] and Diakonov f20fl . For reviews we refer to Refs. |plH?2 



Gradually, since instantons are now successfully identified in lattice ensembles (see the recent 



rapporteur talks |^,p4j]), the instanton vacuum accumulates direct support from lattice QCD. 
However, there are still systematic uncertaincies concerning the parameters of the instanton 
liquid emerging from lattice analyses. 

Without doubt, instantons play an important role in explaining chiral symmetry breaking 
and other empirical facts of hadron structure. However, their contribution to the confinement 
property of non-Abelian gauge theory (as far as this is established on the lattice) turned out 
quantitatively insufficient p5| , p6| , at least at the present level of analytically dealing with the 
instanton liquid model or implementing it in numerical simulation. One step beyond the current 
wisdom, i.e. beyond the neglect of instanton color correlations, has been attempted in Ref. 



27fl. The method proposed there uses the field strength correlator in a dedicated way between 



clusters of topological charge. However, based on the RG smoothing method, the instanton 



interpretation of the emerging clusters in [27] has remained uncertain. Therefore, the same type 
of measurement should be applied in the context of any cooling study devoted to instantons. 

Relating instantons to the confinement issue, alternatively to the QQ force, the property 
of monopole percolation in a multiinstanton system |25| , p8| or the effect of instantons on the 
vacuum structure as summarized in the field strength correlators can be investigated. Earlier 
studies of field strength correlators based on the instanton liquid model can be found in Refs. 
[PP|-|3T|. Here, there are two lines of thought. One P0| , |3"!2|] relies on the one-instanton ap- 



proximation and concentrates on the modification of the single instanton approximation (and 



the instanton solution by itself) due to the interaction with the vacuum medium surrounding 
it. Our approach tries to separate the contributions of one (anti) instanton from contributions 
when the field strength comes from the non-linear superposition of two (anti)instantons. Cor- 
respondingly, different results can be found in Refs. p9|-pT|. 

As far as the one-instanton contribution is concerned, we have nothing to add to our pre- 
vious paper Concerning the two-instanton contributions, we had followed the current 



wisdom and neglected eventual color correlations in the instanton liquid. This is not a bad 
approximation for zero temperature quenched QCD, and we keep this approximation also now. 
However, our results in the previous paper were obtained omitting the Schwinger-line phase 
factors in dealing with the two-instanton contributions. Some (rough) arguments led us to 
expect a negligible systematic error. This was correctly criticized in |[32|| , mainly because of 
some unavoidable consequences. The investigation to be presented here was triggered by this 
discussion and was aimed to improve our previous results avoiding this simplifying assumption. 
This has led us into a rather involved weighted Monte Carlo treatment of the collective coordi- 
nate integration, taking the Schwinger line factors completely into account in the numerically 
evaluated integrand. It turns out that this treatment of the second order contributions partly 
changes our previous conclusions. 

After recalling, for the sake of completeness, the notation for field strength correlators in 
Section II we concentrate on the discussion of the two-instanton contribution in Section III. 



For the first order results we can refer to Ref. |3T| where they were obtained in an analytical 
way. Conclusions will be drawn in Section IV. 

II. THE FIELD STRENGTH CORRELATORS IN THE SEMICLASSICAL 

APPROXIMATION 

The gauge invariant two-point correlators of the non-Abelian field strength are defined as 
V mva {xi -x 2 ) = (0\Ti{G fip ( y x 1 )S(x 1 ,x 2 )G U(T ( y X2)S t (x 1 ,X2)}\0) , (1) 
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where G^p = T a G a ^ p is the field strength tensor and S(xx, X2) is the Schwinger-line phase factor, 
i.e. the parallel transporter necessary to join the field-strength operators at the points x\, X2 
in order to respect gauge invariance. T a denotes the generators of the gauge group SU(N C ). 
The most general form of the correlator compatible with Euclidean 0(4) rotational invariance, 
adequate for zero temperature, is in the notation of Ref. 

£V,iat(z) = ($tJ,u$pa - SprSpv) {^{x 2 ) + ^l^ 2 )) + 

(2) 

, - - r. r. N OVl(X 2 ) 

~\~ [X nXjyO per XnXfjOpy ~\~ X pX a p V X pX \,0 n(j\ 



with x = X\ — X2 and T*(x 2 ) and T>i(x 2 ) representing invariant vacuum structure functions. 
Instead of considering the invariant functions T)(x 2 ) and T>i(x 2 ) entering Eq. (|2|) we shall 
study longitudinal and transverse correlators 

dx 2 

(3) 



V {l (x 2 ) = V(x 2 ) + V,{x 2 ) + x 2 —^, 



V ± (x 2 ) = V(x 2 )+V 1 ( 



x 2 ) 



In particular, if x = (0, 0, 0, x) is Euclidean timelike, we have 



(4) 



i<3 



X 



2\ 



It is easy to demonstrate that T>\ does not contribute to the area law of Wilson loops |§|§. 
In the perturbative regime both invariant functions T> and T>\ behave like 1/x 4 . Only T>\ 
receives a contribution from one-gluon exchange. On very general grounds, the perturbative 
part of T> (which appears at one loop and higher orders) was recently shown to be cancelled 
in the expression for the string tension by higher correlator contributions Here, we shall 
not discuss the perturbative contributions in more detail. Instead we will concentrate on the 
contribution from instantons as exclusive semiclassical configurations representing the non- 
perturbative part of the correlators. 



For comparison we refer to the lattice data published by the Pisa group for the two-point 
field-strength correlators in pure SU(3) gauge theory at T = ||. Since these data include 
measurements for several bare lattice couplings /3, the cooling method seems to be compatible 
with the correct scaling behavior at high (3 of the structure functions. After applying the cooling 
method the non-perturbative parts of the respective structure function have been fitted by an 
exponential function (with the same correlation length) over all distances between 0.1 and 
lfm. This contribution to T> is considerably bigger compared to T>\. This observation points 
towards a dominance of ( ant i) self dual field entering the correlator (0). In the one-instanton 
approximation T>% is exactly vanishing. The structure functions extracted from the cooled 
Monte Carlo data still exhibit a perturbative tail ~ x~ A . In this case, however, the relative 
size of this contribution to T> and T>\ is not understood, in view of the above considerations. 

In the lowest order (in g 2 ) semiclassical approximation for the field-strength correlator 
the Gaussian integral over quantum fluctuations above a single classical field configuration is 
separated from the gauge invariant product of field strengths to be evaluated for this background 
field. What then remains are zero-mode integrations over the appropriate set of collective 
coordinates (summarized as T which characterizes the classical field) with a density function 
M(r), which results from the integration over non-zero mode fluctuations : 

V mu(T {x) = ^ J dTM(T) ■ Tr {G^xu T)S(x 1 ,x 2 ; T)G ua (x 2 ; T)S\x 1 ,x 2 ; T)} , (5) 

where G flu (x;T) is the field strength tensor and S(xi,x 2 ;T) the Schwinger line corresponding 
to configurations A^x, T). To be more specific, we imagine a model of the vacuum state that 



is semi-classically represented by superpositions of N instantons and N anti-instantons IT? 



N N 

A,(x, r) = £ Ap(x- 7i ) + £ A^x; 7j ) • (6) 

i=l j=l 

The ji ( 7j ) denote the collective coordinates of the i-th instanton (j-th anti-instanton), which 
include the positions z i: the group space orientations and the sizes pj. The integration 
measure in Eq. @ is then expressed by 

N N 

dT = Y\_ dji Y\_ i = d Ziduiidpi , d^j = d 4 Zjdujjdpj . 

i=l 3=1 



For practical use we will consider here only instantons and antiinstantons of fixed size p. This 
corresponds to the instanton liquid model invented in |T{| with a delta-like size distribution. 
A more realistic p-distribution with a selfconsistent exponential infra-red cutoff (allowing to 
satisfy low-energy theorems) can be obtained from the assumption that (anti-) instantons repel 
each other at short distances |p|j20| f| 

If the instanton liquid is sufficiently dilute we can approximate the functional integral by 
an expansion in powers of the (anti-)instanton densities = N/V ( fU = N/V ). Then it is 
natural first to try to neglect possible correlation effects due to interactions between instantons. 

Strictly speaking, the superposition ansatz (|j) makes sense as an approximate saddle point 
of the action only if the vector potentials A^, A^ decrease fast enough. This happens when 
the singular gauge expression is used for the (anti-) instanton solutions A^, A^, instead of the 
regular gauge form |[L7||. The existence of a systematic expansion in higher order contributions 



to the measure A^(T) has been proven in Ref. [35 



In |3]J the leading term in an expansion of ([!]) in terms of the density has been discussed 
in detail. In that approximation the field strength correlator is given by the sum of instanton 
(/) and antiinstanton (7) contributions 



= n 4 J d 4 z Tr {G^Xi, j)S(xx, x 2 ; j)G va (x 2 ] l)S\xi, x 2 \ 7)} (7) 

+K,7 -> ^4,7) • 

The integration over the global group orientation of the respective solution is trivial in this 
case and can be omitted. The Schwinger line phase factor is a path dependent matrix in the 
fundamental representation, 



S(x 1 ,x 2 ;z) =P exp^z^ dt x^A^xit); z] 



*Recently, another interpretation has been put forward for such a shape as a manifestation of sup- 
pression by a Higgs like monopole condensate [34|. 
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where the vector potential A M = T a A^, in the case at hand belongs to the single instanton 
source localized at z. Advantage has been taken from the fact that this phase factor, evaluated 
for a straight line has a closed expression which can be easily contracted with the field strengths. 

As mentioned before, due to the vanishing of T>i for purely (anti) self dual fields, the single- 
instanton approximation reads equally for the longitudinal and transverse correlators 

V ll (x 2 )=V ± (x 2 ) = ^n 2 nI(-), (9) 
3 p 

where n = + is the total density of instantons plus antiinstantons. The function I(x/p) 
(for a straight line Schwinger phase factor) is normalized by 1(0) = 1. It has been numerically 
computed in ]31] in order to deal with the final integration over the instanton position z and 
is plotted in Fig. 1. 



III. REVISITING THE SECOND ORDER INSTANTON DENSITY 
CONTRIBUTIONS TO THE FIELD STRENGTH CORRELATORS 

In this section we shall present the numerical integration giving the next order term in an 
expansion in terms of the instanton density. Now we have to consider the field strength from the 
nonlinear superposition of two different solutions A and B, where both A and B can represent 
an instanton or antiinstanton, 

Gfj, p (A, B) = G W (A) + G W {B) + AG W (A, B) , 

(10) 

AG W (A, B) = -i{[A^ B p ] + [B„ A p }} . 

This field strength is plugged into an integral over the (factorized) two-source weight function 
(with fixed sizes p\= pi = p) 

' D fl va {xi,x 2 ) = \ n i A) n i B) I d 4z i I d*z 2 ( dux I du 2 x 

1 A,B=I,I J J J J 

(11) 

Tr 71), B(x l ,^f 2 ))S(xi,x 2 ; 7i, 72) G ua (A(x 2 , 71), B(x 2 , ^S^Xx^x^ 71, 72)} . 



For the purpose of our calculation, the points x\ and x 2 are sitting on the Euclidean time axis 
while the two (anti)instantons are arbitrarily located in 4-d space-time. One has to keep in 
mind, that those formal contributions to Eq. (|TTD, where both the factors G^A, B) at x\ and 
x 2 would receive contributions from only one and the same G^A) (or G^B)) do not really 
occur. They are already taken into account in the first order contribution (|7|). 

With the notation y 1 = x — Z\, y 2 = x — z 2 the two instanton vector potentials, both in the 
singular gauge, read as follows 

2 p 2 

Ki x \ z i) = Va^yiuf{yi,p), B a (x-z 2 ) = uj aa *ri a/ ^ y 2u f(y 2 ,p) , f{y,p) = 2 2 — ^ , (12) 

y \y + p ) 

where u aa i is included to take the relative color orientation inside the pair into account. f} ai iv 



and rj afJLV are the 't Hooft tensors ||T6||. For an antiinstanton instead of an instanton replace 

77 — 5- 77. 

The rotation of instantons in color space (the case of SU(2) group) 

a; - uj ab Al 

can be represented by the related rotation in 3 — d coordinate space 

u ab A^x ,x)=u^A^x Q ,ux), 

u ab A\(x ,x) = A%(x ,ujx) . 

This property will be very important in the further average over the relative color orientation 
of two instantons. 

First, let us consider the correlator T>u(x 2 ). According to Eqs. (10>|lT), the evaluation of 



(fD naively gives rise to 16 terms. After averaging over the relative color and 3 — d space 
orientations of two instantons only six of these terms give a nonzero contribution. Let us 
describe the contributions in some detail. 

1. Two terms, for which both the Gn(A(xi, 71), B(xi, 72)) and the G^(A(x 2 , 71), B(x 2 , 72)) 
are represented by one and the same G^A) (or G^B)) do not really occur (see above). 



2. If in the place of Ga{A(xi, 71), B(xi, 72)) the term G^(A) is taken (or G^(B)) and in the 
place of G±i{A{x2, 71), -6(2:2, 72)) the term Ga{B) (or G 4 j(A), respectively) the resulting 
contributions (2 terms) are identically zero. This can be easily seen. Under color rotation 
of the second instanton B (or A, respectively) (being equivalent to a 3 — d rotation) 
the term in place of G±i{A{x2, 71), B(x2, 72)) becomes rotated, while that in place of 
G ii {A{x 1 , 71), 72)) does not. The rotation of the coordinates of the second instanton 
in the phase factor S(x±, x 2 ; 7172) etc. can be skipped, because there is an integration 
over these coordinates. As a result we have to average over the 3 — d vector G,u(B) (or 
Gn(A)) orientation leading to a vanishing contribution. 

3. If in the place of Ga{A{xi, 71), B(xi, 72)) the term Ga{A) is inserted and in the place of 
Gn{A(x2, 71), B(x2, 72)) the term —i[A±,Bi\, the contribution is zero by the same reason 
as in the previous consideration. There are three other vanishing contributions to be 
obtained from the described one by the replacement A <-> B, x 1 <-> x 2 (4 terms). 

4. If in the place of Ga{A{xi, 71), 72)) the term G 4i (A) is inserted and in the place of 
Ga{A(x2, 71), 5(^2, 72)) the term — i[f? 4 ,Aj] the corresponding contribution is nonzero. 
It is independent of the color orientation of the second instanton B (equivalent to its 
rotation in 3 — d space) as far as the rotation of the second instanton in the phase factor 
S(xi, £2; 7172) etc. is integrated out. There are three other analogous terms with nonzero 
contribution which can be obtained from the described one by the replacements A <-> B, 
x\ <-> x 2 (4 terms). 

5. If in the place of G^{A{xi, 71), B(x±, 72)) the term — i[A 4 ,Sj] (or —i[B 4) Aj\) is inserted 
and in the place of Ga{A{x2, 71), B(x 2 , 72)) the term — i[_B 4 ,Aj] (or — iL4 4 ,-Bj], respec- 
tively), the contribution is vanishing by the same reason as discussed in points 2, 3 above 
(2 terms). 

6. Finally, if in the place of G 4i {A{x 1 , 71), B(x±, 72)) the term — i[A 4 ,Sj] (or — i[_B 4 ,Aj]) is 
inserted and in the place of Ga{A{x2, 71), B(x 2 , 72)) the term — i[A 4 ,Sj] (or — i[S 4 ,Aj], 
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respectively), the contribution is nonvanishing and independent of color rotation of the 
second instanton B (equivalent to its rotation in 3 — d space) by the same reason as in 
point 4 above (2 terms). 

Now let us consider the correlator T>j_(x 2 ). The corresponding expression in (Bl) naively contains 
16 terms upon insertion of ( |T0] ) into Eq. (|TTJ) . After averaging over the relative color and 3 — d 
space orientations of two instantons only four of them give nonzero contribution. 

1. Two terms, for which both the Gij(A(xx, 71), B(xx, j 2 )) and the Gij(A(x 2 , 71), B(x 2 , 72)) 
are represented by one and the same Gij(A) (or Gij(B)) do not really occur (see above). 

2. If in the place of Gij(A(xi, 71), B(x\, 72)) the term Gij(A) is taken (or G{j(B)) and in 
the place of Gij(A(x 2 , 71), B(x 2 , 72)) the term Gij(B) (or Gij(A), respectively) the con- 
tribution is zero. Under color rotation of the second instanton B (or A, respectively) 
which is equivalent to its rotation in 3 — d space Gij(A(x 2 , 71), B(x 2 , 72)) is rotated, while 
Gij (A(xi, 71), B(xi, 72)) is not. The rotation of the coordinates of the second instanton 
in the phase factor S(xi, x 2 ; 7172) etc. can be skipped, because of the integration over 
these coordinates. As a result we have the average over the orientation of a 3 — d vector 
(an antisymmetric tensor is equivalent to a 3 — d vector) Gij(B) (or Gij(A)), which gives 
zero (2 terms). 

3. If in the place of Gij(A(xi, 71), B(x\, 72)) the term Gij(A) is inserted and in the place of 
Gij(A(x 2 , 71), B(x 2 , 72)) the term —i[Ai,Bj] the contribution is zero by the same reason 
as in point 2. There are three other terms with zero contribution which can be obtained 
from the described one by the replacements A «-> B, x\ «-» x 2 (4 terms). 

4. If in the place of Gij{A(x\, 71), B(x\, 72)) the term Gij(A) is inserted and in the place of 
Gij(A(x 2 , 71), B(x 2 , 72)) the term —i[Bi,Aj] the contribution is zero by the same reason 
as in point 2. There are three other terms with zero contribution which can be obtained 
from the described one by the replacements A <-» B, x\ <-> x 2 (4 terms). 
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5. If in the place of G i j{A(xi 1 71), B(x\, 72)) the term — i{[Ai, Bj] + [Bi, Aj]} (at 27) is inserted 
and in the place of Gij(A(x2, 71), B(x 2 , 72)) the same term — Bj] + [Bi, Aj}} (at x 2 ) 
the contribution is non- vanishing (4 terms). 

With the second order contributions in terms of the instanton density included and with an 
additional (unknown) perturbative short-range x~ 4 -contribution the longitudinal and transverse 
correlators can be directly expressed as follows 

^ii(- 2 ) = ^ 2 ^) + ^W/|i(^) + ^, 

(13) 

where n is the total density of instantons and antiinstantons and p is their size, a\\ and a± 
represent the coefficients of the respective perturbative contributions. 

The functions I(x/p), I\\(x/p) and I±(x/ p) have been obtained by numerical integration and 
are all normalized to 1 at x — 0. The numerical factors come from the transition to the SU (3) 
gauge group into which the SU (2)-instantons are embedded and can be calculated analytically 
by the consideration of correlators at x — 0. There is a twofold nontrivial integration in I(x/p), 
and a fivefold integration is implied by the expressions for I\\(x/ p) and I±(x/p). In the last case 
the integrations have been carried out by means of a Monte Carlo importance sampling method. 
The two points x\ and x 2 of the correlator, with x = \xi — x 2 \, are located on the Euclidean time 
axis, at times ±x/2. The two-instanton (or instanton-antiinstanton) contribution is obtained 
by integrating over 77, t 1 , r 2 , t 2 and 9 in a sequential manner, where the positions of the two 
instanton centers are x[ = (0, 0, 77, t\) and x\ = (r 2 sin 9, 0, r 2 cos 9, t 2 ). First, variables (77, ti) 
were generated in a box of size (20 p, 40 p) with a distribution proportional to r\ /(in) /(112), 
where f(x) = | ( x 2 P +p 2) is the profile function of an instanton. Here, in and 272 are the A — d 
distances of the first instanton x[ from 27 and £2, respectively. For this sampling of (77, ti) the 
acceptance varied from 1/1000 to 1/100 for < x < 5 p. About 7000 events were accepted. 
Second, for each accepted event new variables (r 2 , t 2 ) were generated in a similar way, and about 
100 events were accepted per each x[ = (0, 0,77,^1). Finally, for each of the events accepted 
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so far, an angle 9 of relative orientation in 3 — d space has been randomly selected 100 times, 
according to a flat measure in cos 9. The obtained accuracy of the Monte Carlo integration 
was about 1 %. The convergence of integration has been verified by doubling the integration 
box for (r 1 ,ti) and (r 2 ,t 2 ) to (40 p, 80 p). 

The functions Ldx/p) and I±(x/p) are plotted together with I(x/p) obtained in Ref. |31 



in Fig. 1. The longitudinal and transverse correlators from Eqs. (|T3|) have been jointly fitted 
to the lattice data of Ref. |J within a distance range from 0.4 fm to lfm. For the parameters 
n, p, a\\ and a± we have found 

n = (1 ±0.1) far 4 , p= (0.42 ±0.01) fm, a,, = 0.46 ± 0.02 , a ± = 0.76 ± 0.06 

with an xV-^dof = 56.1/(30 — 4). The corresponding curves for the correlations functions T>\\ 
and V± are drawn in Figs. 2 and 3, respectively, together with the lattice data. In order to 
give an impression how large the contribution of the different terms in Eq. (|l^) are, we have 
plotted them separately. 

IV. CONCLUSIONS AND DISCUSSION 

We have considered the two-point correlators of gluon field strengths in the uncorrelated 
instanton liquid model up to the second order in the instanton density. The correlators have 
two parts: a nonperturbative one which is considered exclusively due to instantons and a per- 
turbative one due to gluon exchange fluctuations in the vacuum. We have fitted the resulting 
expression directly to the lattice data within a restricted distance range. For this range the 
achieved quality of our fit is comparable with those of 0, where purely exponential terms 
together with a; _4 -contributions were fitted to the lattice data after a separation into T> and 
T>\. The value for the instanton density comes near to the value expected from phenomeno- 
logical applications, although our present analysis describes quenched lattice data. There is a 



recent analysis due to A. Hasenfratz |36| which, on the basis of the two-point correlator of the 
topological density, gives estimates for the instanton density, the fraction of two-instanton and 
instanton-antiinstanton molecules, all for quenched SU(3) theory and full QCD. The instanton 
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density was there found to be about 1 fm -4 for pure SU (3), too, however the instanton size was 
estimated as about 0.3 fm. Thus, the mean instanton size obtained here is more in accordance 
with UKQCD |37[]. The packing fraction of the instanton gas is estimated to np 4 ~ .03 which 



corresponds to a reasonable diluteness. 

In our fit for the range 0.4 fm < x < 1 fm the instanton and perturbative contributions 
are of a comparable order of magnitude. We obtain the relation V± > Du now both for 
the perturbative and the instanton contributions. This means that both contributions to the 



invariant function V>\ turn out to be positive. This is contrary to our previous result in pi 
where we concluded that T>\ < is inavoidable within the instanton liquid description. The 
neglect of the Schwinger-line phase factor in the two-instanton (instanton-antiinstanton field) 
is to be blamed for this wrong conclusion. 

The second order density terms are the final ones for transverse correlators, i.e. the 
0{n z ) ) 0{n i ) 1 - ■ ■ terms would be vanishing. For the longitudinal correlators there exist also 
a non- vanishing 0(n 3 ) contribution, but no 0(n 4 ) ones. The calculation of the longitudinal 
correlators up to this order require a nine- dimensional nontrivial integration over the three- 
instanton degrees of freedom. We hope to come back to this integration in future. 
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FIGURES 




FIG. 1. The numerically computed functions I(x/p) for the first order instanton contribution (from 
[31]) and I^(x/p), I±(x/p) for the second order contribution. 
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FIG. 2. Transverse field strength correlator as described by the instanton gas model (continuous 
line) and fitted to the lattice data points of Ref. [9] The dashed lines show separately the first and 
second order instanton density contributions as well as the perturbative x~ 4 contribution. 
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